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Both simple and sophisticated models are frequently used 
in an attempt to understand how real nuclei breakup when 
subjected to large excitation energies, a process known as 
nuclear multifragmentation. Many of these models assume 
equilibrium thermodynamics and produce results often inter- 
preted as evidence of a phase transition. This work examines 
one class of models and employs standard thermodynamical 
procedure to explore the possible existence and nature of a 
phase transition. The role of various terms, e.g. Coulomb 
and surface energy, is discussed. 



I. INTRODUCTION AND OVERVIEW 

Many models |jj - M have been proposed the describe 
the breakup of a large nucleus subjected to excitation 
energies greater than a few MeV per nucleon, a process 
known as multifragmentation. Experimentally, the sig- 
nature of multifragmentation is the production of a wide 
range of nuclear reaction products, particularly interme- 
diate mass fragments (IMFs), 3 < Z < 30. On the basis 
of inclusive data, it was proposed [^0| - [Q that these 
fragments were produced in analogy to a liquid-to-gas 
phase transition occurring in a nucleus. A recent exper- 
iment that permitted the total charge reconstruction of 
each event studied multifragmentation resulting from the 
breakup of gold nuclei as a function of the excitation en- 
ergy deposited Q - (lTj . The statistical aspects of these 
data have provided strong evidence that multifragmenta- 
tion is indeed related to a phase transition occurring in a 
finite system. Whether the production of IMFs in such 
collisions is due to a phase transition, and if so, what 
type, is still an issue of much debate fl8[| . 

One class of models developed to explore the fate of a 
nucleus as a function of excitation energy is based on the 
phenomenological description of the free energy, F(V, T), 
of the breakup state, where T is the common temperature 
of all nucleons and nuclei within the breakup volume V. 
These nuclei are considered to be at normal nuclear den- 
sity and interact only via the Coulomb force. The distri- 
bution of nuclear fragments prior to any secondary decay 
can then be calculated as a first step in the disassembly 
of the excited initial system. To compare with data, de- 
excitation of fragments and expansion of the system due 



to the Coulomb repulsion between the fragments must be 
accounted for in the model. However, if the thermody- 
namics of the model is of interest, as is the case in this 
work, then only the behavior of thermodynamic variables 
need be examined, e.g. free energy, entropy, specific heat, 
pressure, isothermal compressibility. Thus, no fragment 
distributions need be explicitly calculated and therefore 
no fragment distributions are analyzed in this work. 

Here, several variations of a previously discussed model 
[pi are explored. A canonical ensemble approach is used 
to investigate the thermodynamics of the system where 
the free energy, F, is written as a function of the tempera- 
ture, T, and the volume, V . Calculations are restricted to 
a system which contains 162 constituent, since this is rep- 
resentative of the size of the system studied in . Con- 
tributions to F, e.g. the surface free energy, the Coulomb 
energy, are examined by turning them off or altering the 
form of the contribution in question. In this way insight 
can be gained as to how the important features of the 
thermodynamics such as specific heat, isothermal com- 
pressibility, etc. depend on the parameterization of the 
free energy. 

This paper is organized as follows. In section II the 
details of the models are presented. Three versions of 
a standard statistical multifragmentation model are ex- 
amined as well as a well-known mean field model whose 
results are used for comparison. In section III a descrip- 
tion of the analysis and the results of that analysis are 
presented. Section IV discusses the standard interpreta- 
tions of the models and analysis. Finally, a brief discus- 
sion of the questions raised by this work concludes this 
paper. In general, the notation of references 0] and ||] 
are followed. 



II. DETAILS OF THE MODELS 

This work follows directly the efforts presented in ref- 
erence H in which the canonical partition function was 
examined as a function of temperature in a fixed volume 
system for evidence of a phase transition. In that work, 
evidence for a first order phase transition was found. In 
the present work, the volume (average density) of the 
system is permitted to vary. It shall be seen that the 
nature of the phase transition depends on the volume of 
the system. The work of || is also extended by exam- 
ining the effects of the Coulomb force on the system and 
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the effects of the choice of surface energy parameteriza- 
tions. The units used for the nuclear models are: energy 
and free energy in MeV/nucleon, temperature in MeV, 
volume in fm 3 , pressure in MeV/fm 3 and so on. 
A general description of each system follows. 

A. VI: Full statistical description of an excited 
nucleus 

Calculations begin by considering the free energy of a 
nuclear fragment. It is assumed that the free energy 
of a nuclear fragment of mass A and charge Z, for A > 1 
is given by: 

Fa.z = F B KZ + F%? + Fl z + E% z . (1) 

The terms in eq. (|l|) refer to the bulk, symmetry, surface 
and Coulomb contributions to the free energy of a nuclear 
fragment. The forms of these terms are given Q: 

Fi z = {-W -T*/e )A, (2) 
FZz=l{A-2Zf/A, (3) 

/rp2 _rp2 \ 5 / 4 

E C A .z = \e 2 Z 2 {l (1 + «)-V3)/i2 ( 5 ) 
5 

In eq. (||) the constants are taken as Wq = 16 MeV and 
e Q = 16 MeV. In eq. © 7 = 25 MeV. In eq. (J§) 
00 = 18 MeV and T c = 16 MeV, following reference j|. 
The contribution from the Coulomb term is estimated 
via a Wigner-Seitz approximation as in reference Q|. 

The re-term is related to the volume of the system 
through 

1 + k = V/V . (6) 

This simplified model presented here differs from the 
standard version jjj in that there is only one parame- 
ter relating the volume excluded by the constituents, Vo, 
to the total volume of the system, V, and to the free 
volume, Vf. Here the free volume is the difference be- 
tween the total volume V and the sum of the volume of 
the fragments, assumed to be at normal nuclear density, 
and is the volume available for the translational motion 
of the fragments. 

In the standard version of the model the free volume is 
given by Vf = "xYa , where x is parameterized to increase 
with fragment multiplicity such that it varies between 
0.2 and 2; the parameter re is fixed, usually at re = 2. 
For simplicity, here it is assumed that re = x so that 



specifying Vf determines the value of re in eq. (JsJ) . See 
reference jlj for details of re and x m t ne standard version. 

For this work then, the total volume of the system is 
then given by: 

V = Vo + V f . (7) 

Two things become obvious from eq. (Q); first, with this 
form of V the free energy of the system varies with Vf 
since Vo is a constant. Second, the loss of free volume in 
the closest packing of spherical clusters is ignored. The 
issue of whether spherical nuclei can actually be placed 
in a total volume, Vo, given a free volume, Vf, is not 
addressed. Undoubtedly, there will be situations where 
it is not possible for the total volume to accommodate 
all of the nuclear clusters. The purpose here, however, 
is to explore the thermodynamics and self-consistency of 
the model and not physical consistency. 

Finally Ra,z is the radius of the fragment in question 
and is determined by 

Ra,z ^roA 1 / 3 , (8) 

with ro = 1.17 fm. The version of the model presented 
above will be termed VI. 

From this point the intrinsic partition function of a 
fragment of A, Z at temperature T and volume V can 
be determined ||] as follows: 

z A , z = exp(-F A ^/T). (9) 

Using a technique developed in reference [jl9| and used 
on a simplified version of this model B, the canonical 
partition function can be built via a recursion relation: 

1 P 

Z p = - ^ Au! A Zp_ A , (10) 

p A=l 

starting from Zo — 1. Here for calculational simplicity 
the approximation has been made that for each and every 
fragment with A > 1, A/Z — 2.5 which represents an 
average mass to charge ratio for fragments. The lo a term 
is 

, a = ^(2™T) 3 / 2 # 2 z,, z , (11) 

where the terms the the left of the fragment partition 
function, z A ,z, account for the translational free energy 
contribution, F tr . It is now straightforward to calculate 
the partition function of the system for a given T, Vf, Aq 
and Zo- The free energy of the system of j?-particles is 
then determined as usual 

F = -Tlog(Z p )+E^V), (12) 

where the last term is the usual Coulomb contribution of 
a uniformly charged sphere: 

W) = P^ (13) 
with R = (3U/47T) 1 / 3 . 
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1. Comparison of VI to the full version of the model 

The model and calculations described above were com- 
pared to the full, or unmodified version of the model of- 
ten cited in the literature, see for example references [Q, 
§0j - ||||. See Figure 1. In Figure 1 results from the 
full version of this model are shown for the mean frag- 
ment distribution calculated at a given input excitation 
energy. To generate event-by-event distributions Poisso- 
nian fluctuations about the mean are introduced, after 
which, temperature is adjusted to ensure energy conser- 
vation. To more fully recover the standard version of 
the model most often used, higher order corrections were 
introduced just as in the full version of the model; e.g. 
eo in cq. (|^) was made dependent on the fragment mass 
A, for light clusters, A < 4, the empirical masses and 
binding energies, radii and spin degeneracy factors of the 
ground state were used, the total volume was held con- 
stant at 3Vq and the free volume was set to depend on 
the input excitation energy. Finally, energy was explic- 
itly conserved; an input excitation energy was given and 
a temperature was determined such that total energy was 
conserved. The explicit conservation of energy produced 
results that were essentially the same as those resulting 
from the unconstrained canonical ensemble. 
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FIG. 1. Comparison of the full model (open circles) before 
Poissonian smearing and the calculations of this paper (full 
circles) for the caloric curve (a) and the average fragment 
multiplicities (b). 

Figure la shows the caloric curve from the full version 
of this model for a system with 100 nucleons (60 neutrons 



and 40 protons) compared to the same size system used in 
calculations with a modified version of VI. The general 
trend of the modified VI reproduced the average behavior 
of the full model, though there is not a perfect agreement. 
This is to be expected. While this modified version of 
VI is closer to the model, there are still some differences, 
e.g. the charge of A > 4 fragment is treated in only an 
average fashion in VI. The reproduction of the general 
trends indicates that VI captures the essence of the full 
model. Figure lb shows the fragment multiplicity, before 
any secondary decay, from both models. Again there is 
general agreement between the two. 

The break observed in the caloric curve shown in Fig- 
ure 1 is well known in the full model, see, for example, 
Figure 4 in ref . [^(J and Figure 1 1 in ref . [^|| . The break 
is due to the initial guess of the system's multiplicity 
which is in turn used to guess the system's free volume. 
For low energies the multiplicity is chosen to be 1, 2 or 3 
(Figure lb shows that the initial guess of the multiplicity 
is consistent with the final state multiplicity), while at 
higher energies the multiplicity depends smoothly on a 
function of the input excitation energy jl[] . In some sys- 
tems, e.g. a system of 100 nucleons, there is a jump in 
multiplicity at the transition from the low energy com- 
putations and the high energy computations which gives 
rise to a jump in the final state multiplicity, Figure lb, 
and a break in the caloric curve, Figure la. When the 
simplified model used in this work is given the same vol- 
ume dependence as the full model, the results of the full 
model are reproduced. 

Energy conservation is explored in Figure 2. Here the 
unmodified version of VI was used with a system of 162 
particles and energy was explicitly conserved as outlined 
above. In order to examine the change in energy between 
the initial and final state of each term contributing to the 
system's total energy, the temperature of the initial state 
of the system has been calculated corresponding to the 
input E* . The assumption that the initial nuclear state 
is in thermal equilibrium prior to its deexcitation to the 
final state has no bearing on the thermodynamics of the 
final state and is done only for purposes of the above- 
mentioned calculation. The initial state used for this 
calculation was the system of 162 nucleons at excitation 
energy, E* , at normal nuclear density, po and at a tem- 
perature that conserves energy when the total energy is 
determined using eqns (||)-(||), JT^ ) and: 

E = F - T {§) V <"» 

The total energy of the initial state is shown in Figure 
2 as well as it's various components and a caloric curve. 
The final state of the system was computed with the 
same E* but was held at a third normal density, po/3, 
and allowed to fragment in the manner outlined above. 
The caloric curve produced for the final state via this 
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explicit conservation of energy calculation is identical to 
the caloric curve produced via the calculations without 
an explicit conservation of energy. Figure 2 shows that 
the temperature in the final state is lower than that in the 
initial state. Further inspection indicates that while the 
Coulomb energy is reduced by creating smaller charged 
nuclei, the energy required to create the additional sur- 
face area more is more than offsetting. Thus, the tem- 
perature must decrease. 




FIG. 2. Energy conservation in VI, overall, component by 
component (total, bulk, surface, translational, symmetry and 
Coulomb energies respectively) and the caloric curves. Full 
circles represent the initial state, open circles represent the 
final state and open squares represent the difference between 
between the two. 



C. V3: Description of an excited charge-free nucleus 
with a temperature independent surface 

In this version, the Coulomb force is suppressed and 
the temperature dependent surface term in eq. (|j) is 
made independent of temperature: 



(3 A 2/3 . 



(17) 



The free energy of a fragment and the entire system are 
still given by eqns (15) and (filf). Every other aspect of 
the model is the same as in VI and V2. This version of 
the model will be termed V3. 



D. The van der Waals fluid 

The free energy of the van der Waals fluid is deter- 
mined in the standard textbook fashion. Starting from 
the free energy p3fl 

F/A = -t {In [n Q {V - Ab) /A] + 1} - Aa/V, (18) 

with t = kf,T, hq = (Mt/2irh 2 ^ 3 ^ and a and b the 
usual van der Waals constants. In this work, eq. (|l^) is 
computed in terms of the density, p — A/V, so that the 
number of constituents, A, is not a factor. For the van 
der Waals' constants a and b, values were used for helium 
so that T c ~ 4.5 x 1CT 4 eV/A. This also suggested the 
value of M in tiq. Finally, eq. ( |l8| ) was evaluated in 
terms of T/T c and V/V c . The van der Waals fluid model 
is well defined ad free of internal inconsistencies. It will 
be used as a benchmark for the analysis presented in this 
paper. Units for the van der Waals fluid results will be in 
cV/A for energy and free energy, eV/A x V c for pressure 
and so on. 



III. DISCUSSION OF CALCULATIONS 



B. V2: Description of an excited charge-free nucleus 

The general ideas of VI are followed but with the 
Coulomb energy of the system set to zero. The free en- 
ergy shown in eq. ([!]) then becomes: 

Fa — Fa + Fl (15) 

And the total free energy of the system is given by 

F = -T\og{Z p ). (16) 

Every other aspect of the model is the same as in VI. 
This version of the model will be termed V2. 



Calculations were performed for each version of the 
model to determine F(T, Vf,A n , Z ) for A = 162, 
Aq/Zq = 2.5 and over a range in temperature, 1 MeV 
< T < 14 MeV, and volume of 2 x 1CT 8 < (V f /V ) < 
2 x 10 8 . Once the general vicinity of the critical point was 
identified, a smaller range in (T, Vf) was used for more 
detailed calculations. For the van der Waals fluid the 
range was smaller and in terms of reduced temperature 
and volume: 0.1 < T/T c < 2.0 and 0.34 < V/V c < 2.0. 
Figure 3 shows the behavior of the free energy over the 
ranges of temperature and volume used in the calcula- 
tions. 

After the value of the free energy was calculated, it 
was simple to determine other thermodynamic quanti- 
ties. Holding the volume fixed the entropy is given by 
the usual relation: 
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(24) 




where (Vf) is the average value of Vf over the AVf in- 
terval. With this information it is possible to determine 
if there is a phase transition in a model such as this and, 
if present, the nature of that phase transition. The fol- 
lowing section addresses this question. 

IV. RESULTS OF CALCULATIONS 

In this section, each of the axes of the standard phase 
diagram, T, V and P, will in turn be held fixed. The 
behavior of other quantities will be examined in order to 
understand the behavior of each system. The van der 
Waals fluid will serve as a guide for the interpretation 
of the analysis and also as a benchmark to illustrate the 
accuracy of this analysis. 

A. Isotherms 



FIG. 3. Free energy surfaces as a function of temperature 
and volume for (a) VI, (b) V2, (c) V3 and (d) the van der 
Waals fluid. 
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9Tj Vf 



(19) 



In the case of this work differences were used instead of 
derivatives due to the numerical nature of the calculation, 
thus eq. (Il9|) becomes 



(at) 



(20) 



v f 



Similarly the specific heat at constant volume was deter- 
mined via 



(T) 



[at, 



(21) 



where (T) is the average value of T over the AT interval. 
Using the entropy and the free energy, the total energy 
can be determined from: 



E = F + TS. 



(22) 



In these calculations it was possible to hold either the 
temperature or the volume constant. The pressure was 
then found by holding the temperature fixed 



(23) 



Taking another derivative then gave the isothermal com- 
pressibility 



Determination of the critical point, coexistence and 
spinodal curves is discussed below. The variation of the 
free energy as a function of density shows the same gen- 
eral features for these systems. 




-23.5 
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FIG. 4. Free energy isotherms as a function of reduced den- 
sity for (a) VI, (b) V2, (c) V3 and (d) the van der Waals fluid. 
Open squares show the spinodal, full circles show the bound- 
ary of the coexistence region. Dotted curves are sub-critical. 
Solid curves are critical. Dashed curves are super-critical. 
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Figure 4 shows the behavior of the free energy 
isotherms for each system as function of reduced den- 
sity. Each plot in Figure 3 shows three isotherms. The 
isotherms are sub-critical (T = 0.95T C ), critical (T = 
T c ) and super-critical (T — 1.05T C ). Determination of 
the critical point is discussed below. Also shown are 
the approximate location of the coexistence and spin- 
odal curves. Determination of these curves is also dis- 
cussed below. The behavior of the free energy for these 
isotherms for all systems is more or less the same. As 
the reduced density increases, the free energy increases. 
At some mid range in reduced density the slope of the 
increase in free energy changes. At a greater reduced den- 
sity the slope of the increase in free energy changes again. 
This is most clearly demonstrated by the van der Waals 
fluid system. See Figure 4d. However, the behavior is 
present in all the models. This behavior, while appearing 
modest in these plots, will be seen to be the cause of the 
critical-like behavior exhibited by these models. 




FIG. 5. Isotherms or pressure as a function of reduced den- 
sity for (a) VI, (b) V2, (c) V3 and (d) the van der Waals 
fluid. Open squares show the spinodal, full circles show the 
coexistence curve. The solid horizontal line is an example 
of a Maxwell equal area construction to determine the co- 
existence region and find the critical point. Dotted curves 
are sub-critical. Solid curves are critical. Dashed curves are 
super-critical. 

The pressure was calculated from the free energy 
isotherms via eq. (|2^). Figure 5 shows the results for 
each model. The determination of the location of the 
critical point and the coexistence and spinodal curves is 
based on the phase diagram of pressure, temperature and 



reduced density. By searching for inflections points along 
the pressure versus reduced density isotherms the spin- 
odal curve was determined. The isotherm immediately 
following, as the temperature of each isotherm increases, 
the last isotherm with two inflection points was labeled 
the critical isotherm. 

Another method to determine the location of the criti- 
cal point began with the isothermal compressibility which 
was calculated with eq. (p4]). Isotherms of kt versus re- 
duced density were inspected. All sub-critical isotherms 
showed at least one negative value of kt- The first 
isotherm, as a function of increasing temperature, which 
showed only positive values of kt was labeled the critical 
isotherm. Both procedures yielded the same results, as 
they are essentially identical. See Table I. 

The coexistence curve was determined by making a 
Maxwell equal area construction for each isotherm. See 
Figure 4 for examples of the Maxwell construction. 

On first glance at Table 1 several noteworthy features 
stand out: (1) none of the values of T c determined for 
the models is the same as the value of the parameter T c 
specified in the surface term in eq. (Q); (2) the critical 
densities for each model is close to unity; and (3) the 
critical temperature of VI, the model which includes the 
Coulomb force, is larger than the critical temperature 
of V2, the model with no Coulomb force. The error of 
the analysis of the van der Waals fluid was on the order 
of a few percent; T^ dw ^ 1. This illustrates the error 
inherent with this type of analysis. 

It is not surprising that the critical temperature found 
by the analysis of thermodynamical quantities is different 
than the parameter T c used to parameterize the surface 
free energy of infinite nuclear matter. If one considers the 
critical point to be that temperature at which the surface 
free energy vanishes, then this can only be at T c (16 MeV 
in this case). However, the form of the surface term given 
by eq. (^), approximately the form of the macroscopic 
surface free energy of a fluid near its critical point, is not 
an appropriate description of the microscopic surface of 
a droplet (24]. Moreover, eq. (|J) leads to a specific heat 
which approaches negative infinity as T approaches T c . 
In a more fundamental model the critical temperature 
would be an output of the model rather than an input. 
Of the phenomenological models studied here, only the 
van der Waals model is known to be self-consistent. 

The high densities found at the critical point for VI, 
V2, and V3 cannot be realized if one is constrained to 
placing spherical nuclei without overlaps inside of the 
total volume. Again this issue is mentioned but not dealt 
with since the aim of this paper is to explore only the 
thermodynamical predictions of the above models. 

Were this sort of model interpreted physically, the high 
value of the critical density determined here would sug- 
gest that the critical point could never be reached by fi- 
nite nuclear matter as a multiplicity of spherical clusters 
at normal nuclear density could not physically fit into 
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TABLE I. Critical points for the models 



System 


T c 


Pc 


Pc 


c f 


VI 


7.575 ± MeV 


0.94±0.01 x po 


(3.2 ± 0.1)xl0 _3 MeV/fm 3 
(7.3±0.1)xl0~ 4 MeV/fm 3 


0.49 ±0.02 


V2 


6.875 ± MeV 


0.74±0.01 x po 


0.16 ±0.02 


V3 


10.975 ± MeV 


0.85±0.01 x po 


(2.4±0.1)xl0~ 3 MeV/fm 3 
(1.7±0.1)xl0" 4 eV/AxV c 


0.28 ±0.02 


vdW 


1.01 ±0.01 x T c 


1.01±0.01 x p c 


0.37 ±0.02 



the critical volume. If the constraint that the breakup 
volume must be large enough to avoid overlapping vol- 
umes of the final state (spherical) nuclei is added, then 
only first order phase transitions are possible. Of course 
there are several problems with a strictly physical inter- 
pretation of models such as the one presented here, not 
the least of which is the introduction of a volume for the 
system. Actual nuclei excited to high energies in nucleus- 
nucleus collisions do not exist in a box and thus have no 
volume in the sense suggested here. 

Item (3), the rise in T c with the vanishing of the 
Coulomb force, is, on the surface, counterintuitive. Many 
other models of nuclear systems show just the opposite 
behavior ^6|, ^7j. However, those models are funda- 
mentally different than the ones examined in this work. 
Such models begin by describing the free energy or chem- 
ical potential or some equivalent quantity using formulae 
which assume a uniform distribution of material in much 
the same way that the van der Waals fluid assumes a 
uniform density. The model in this work samples a dif- 
ferent part of the final state phase space |^5| . It will be 
seen that the calculation of the Coulomb term via eq. (||) 
gives rise to the counterintuitive rise in the critical tem- 
perature when the Coulomb force is suppressed in the 
model. 

To begin to understand the effect of the Coulomb en- 
ergy on the critical point the same isotherm for models 
VI and V2 was examined. Figures 4a and b show the 
isotherm of T = 7.2 MeV for VI and V2. When the co- 
existence curve is shown for both systems, it is obvious 
that for VI the isotherm is sub-critical while for V2 the 
isotherm is super-critical. Figures 4c and d begin to shed 
light on the cause of this counterintuitive occurrence. 

For the sake of illustration, the Maxwell constructed 
(a line of constant slope through the coexistence region, 
which then leads to a constant value of the pressure 
through the coexistence region that will give "equal ar- 
eas" on a pressure-volume plot, see Figure 5a) path of 
the free energy through the coexistence region is shown 
as a dashed line, barely visible just below the isotherm 
in the coexistence region, in Figure 6a. The Maxwell 
constructed free energy is a straight line through the co- 
existence region. The constant slope of the Maxwell con- 
structed free energy leads to a constant pressure for the 
system in the coexistence region. Because the path of 
Maxwell constructed free energy is very close to the path 




P/p„ P/P„ 



FIG. 6. Free energy isotherms as a function of reduced den- 
sity for (a) VI and (b) V2. For both systems the isotherm 
is for T — 7.2 MeV. Open squares show the spinodal, full 
circles show the coexistence curve. The Maxwell constructed 
free energy is shown as a dashed line through the coexistence 
region. The difference of the computed free energy and the 
Maxwell constructed free energy shows the presence of inflec- 
tions points in (c) VI at the coexistence boundary for this 
sub-critical isotherm and none for (d) V2 for this super-critical 
isotherm. 

of the free energy of VI, it is difficult to see the difference 
in a plot such as shown in Figure 6a. A plot of the differ- 
ence in the calculated free energy of VI and the Maxwell 
constructed free energy shows what gives rise to the van 
der Waals loops in Figure 5a. See Figure 6c. There are 
two inflection points in the curve of the calculated free 
energy of VI, these are shown more clearly as the points 
in Figure 6c where the curve shows an ordinate value of 
zero. These plots for a canonical system's free energy are 
in the same spirit as plots for a microcanonical system's 
entropy [p8j . It is clear from Figures 4b and d that there 
are no similar inflection points in the free energy curve 
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for V2. Therefore the isotherm of T = 7.2 MeV in V2 is 
super-critical while the very same isotherm is sub-critical 
for VI. 




0.7 0.8 0.9 1 0.7 0.8 0.9 



P/P„ P/P„ 

FIG. 7. Components of the free energy isotherms (T = 7.2 
MeV) for VI (left) and V2 (right). From top to bottom the 
components are for translation, bulk or Fermi, surface, sym- 
metry, total Coulomb, uniform background Coulomb contri- 
bution and the clusterization Coulomb contribution. Dotted 
lines show the coexistence boundary for this isotherm for VI 

It is possible to understand what gives rise to the in- 
flection points introduced when going from V2 to VI by 
looking at the contributions to the total free energy of 
each system along the isotherm T — 7.2 MeV. See Fig- 
ure 7. Plotted in this figure for both systems are the 
components of the overall free energy for the isotherm 
in question: translational, bulk, surface, symmetry, to- 
tal Coulomb, background Coulomb and clusterization 
Coulomb free energies. Figure 7 shows that the trans- 
lational and bulk free energies of each system are nearly 
identical. An inspection of eqns ([ll|) and (||) shows that 
these quantities are relatively insensitive to small changes 
in the fragment distribution. On the other hand, the sur- 
face free energy shows an obvious difference in behavior 
between systems. In VI the initial decrease in surface free 
energy as a function of reduced density is slower than in 
V2. An inspection of eq. (ji|) shows that the A 2 / 3 -term 
introduces more sensitivity to the fragment distribution 
than the previously discussed terms. The cause of this 
difference, and all the differences between these two sys- 
tems, is the presence of the Coulomb force in VI and its 
absence in V2. 

The behavior of the Coulomb contribution to the free 



energy of this isotherm is now discussed. The first, and 
simplest, additional term in the free energy due to the 
Coulomb force is the asymmetry term. The change is in 
the asymmetry term is smooth as a function of increasing 
density and will not introduce the inflection points in 
the free energy curve that will change a super-critical 
isotherm to a sub-critical isotherm. 

In this model the total Coulomb contribution to the 
free energy comes from two sources. One which repre- 
sents the background energy due to a uniform distribu- 
tion of charges, eq. (|l3|), and the other due to the energy 
from the clusterization of fragments, eq. (@) jlj. The 
background energy, E§, goes as p 1 / 3 and therefore varies 
smoothly with volume. The clusterization free energy, 
Ejlz-i shown in Figure 7 has a different behavior. It 
decreases as the reduced density increases, the rate of 
decrease is at first nearly constant, then slows and then 
increases rapidly over some small interval in reduced den- 
sity. The combination of these two Coulomb terms in- 
troduces sufficient changes in the overall free energy of 
the system from V2 to VI that inflection points arise and 
thus the critical temperature increase when the Coulomb 
force is added to the system. A smooth or constant ver- 
sion of the Coulomb free energy added to version V2 
should not cause this sort of behavior. It is the varia- 
tion in E^ z that introduced the inflection points and 
increases the critical temperature. In the end, it is the 
behavior of the free energy curve that served to deter- 
mine the location of the critical point and that behavior 
is, at times, counter intuitive. 

Also listed in Table I is the compressibility factor: 

P V 

C f = ^f- (25) 

The text book value for Cf for the van der Waals gas is 
recovered to within error bars. According to the law of 
corresponding states, the value of Cf should be univer- 
sal. For fluid systems this is the case and (C/) ~ 0.292 
p9| . For VI and V2 there no such universal behavior ob- 
served to within error bars, while V2 shows some degree 
of universality. 

B. Isochores 

The volume of the system is held constant and the be- 
havior of various quantities with respect to the system's 
temperature is explored. 

Beginning again with the primary quantity calculated, 
Figure 8 shows the free energy for each system as a func- 
tion of the system's temperature. For models VI, V2 
and V3 the isotherms are for p — po/3, p — p c , and 
p = 0.9995/90- For the van der Waals fluid, the isochores 
shown are for p = 1.25p c , p = p c and p — 0.75p c . See ta- 
ble I for critical density values. Also shown are the values 
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FIG. 8. Free energy isochores for (a) VI, (b) V2, (c) V3 and 
(d) the van der Waals fluid. Full circles show the boundary of 
the coexistence region. Dotted curves are sub-critical. Solid 
curves are critical. Dashed curves are super-critical. 



of the free energy of the systems along their respective 
coexistence curves. 

In reference || the authors show results of a similar 
model similar to V2, for an isochore of approximately 
P — Po/3- In that work there was a kink in the free energy 
curve which was interpreted as evidence for a first order 
phase transition. For the small system used in this work 
the kink is smoothed out negating the efficacy in using 
the kink as evidence towards determining the order of 
a phase transition if one is present. In larger systems, 
systems in ref. [|J were more than 15 times larger than 
the system used in this work, the kink is more evident 
and this procedure may be possible. 

Knowledge of the location of the coexistence curve 
allows for the identification of sub-critical, critical and 
super-critical isochores. For the projections shown in 
Figure 8 all the versions of the nuclear model a sub- 
critical isochore crosses the coexistence region. The criti- 
cal isochore travels along the high value edge of the coex- 
istence curve passing through the critical point. Super- 
critical isochores do not come into contact with the coex- 
istence curve nor do they traverse the coexistence region. 
The behavior of the van der Waals fluid is somewhat dif- 
ferent from the behavior of the nuclear models. Figure 
8d shows that for a van der Waals fluid in this projec- 
tion of the phase diagram all are in the coexistence re- 
gion for low temperatures and cross the coexistence curve 
at higher temperatures. This apparent behavior results 



from projection of the three dimensional phase diagram 
onto a two dimensional plot. In a three dimensional fig- 
ure, the super-critical isochore is seen to travel outside 
the coexistence region, the critical isochore is observed 
to intersect with the critical point and the sub-critical 
isochore is seen to traverse the coexistence region. Sec 
Figure 9. 
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FIG. 9. A three dimensional plot of the free energy iso- 
chores for the van der Waals fluid. Full circles show the 
boundary of the coexistence region. Dotted curves are 
sub-critical. Solid curves are critical. Dashed curves are su- 
per-critical. 

The pressure isochore was obtained by following the 
procedure outlined in eq. (23) and using (p), the aver- 
age density, and (T), the average temperature over the 
AT interval. See Figure 10. In the limit of vanishing 
AT this procedure is valid. Figure 10 shows the same 
three isochores discussed above as well as the location 
of the coexistence curve, seen edge on in this projec- 
tion of the phase diagram. All the figures for the nu- 
clear model show similar behaviors for the critical and 
super-critical isochore. All super-critical isochores fol- 
low a trajectory above the coexistence curve. The criti- 
cal isochore for all versions of the model approaches the 
coexistence curve, follows it, and leaves at the termina- 
tion or critical point. The sub-critical isochores for V2 
and V3 pass through the coexistence curve as the tem- 
perature is increased and the system makes a first order 
transition from a liquid to gaseous state. For version VI 
the behavior is more complicated. The sub-critical iso- 
chore begins on the gaseous side of the coexistence curve, 
crosses the coexistence curve into the liquid region before 
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crossing back over the coexistence curve into the gaseous 
region at higher temperature. The P-T projection of the 
van der Waals fluid looks as expected. See Figure 8d. 



the various isochores. 
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FIG. 10. Isochores of pressure versus temperature for (a) 
VI, (b) V2, (c) V3 and (d) the van der Waals fluid. Full circles 
show the boundary of the coexistence region. Dotted curves 
are sub-critical. Solid curves are critical. Dashed curves are 
super-critical. 



Following the analysis procedure outlined in the previ- 
ous section, the entropy was determined via eq. (20). 
See Figure 11. The same three isochores are plotted 
showing the entropy as a function of temperature. All 



three isochores for the nuclear models show a smooth 
rise as a function of temperature with some region of in- 
creased slope. This behavior is consistent with a contin- 
uous phase transition, if a phase transition were present. 
Along an isochore, a first order transition would be indi- 
cated by a sudden change the behavior of the entropy as 
a function of temperature at one edge of the coexistence 
region. However, due to the small size of the system such 
sharp behavior is smoothed out into smooth curves mak- 
ing it impossible to draw a conclusion about the order 
of a phase transition from plots such as those shown in 
Figure 11. As before, knowing the location of the co- 
existence curves makes identification sub-critical, critical 
and super-critical isochores in the nuclear model trivial. 
The sub-critical isochore traverses the coexistence region, 
the critical isochore passes through the critical point and 
the super-critical isochore avoids the coexistence region. 
Also as before the van der Waals system shows a dif- 
ferent behavior, and an added dimension to the plot in 
Figure lid must be made to understand the behavior of 




FIG. 11. Isochores of the entropy versus temperature for 
(a) VI, (b) V2, (c) V3 and (d) the van der Waals fluid. Full 
circles show the boundary of the coexistence region. Dot- 
ted curves are sub-critical. Solid curves are critical. Dashed 
curves are super-critical. 

The specific heat at a constant volume, Cy, for each 
system is shown in Figure 12, again for the same three 
isochores. All curves show a peak in that could be due 
to a smoothed out discontinuity (first order phase tran- 
sition), the remnants of a power law divergence at the 
critical point (continuous phase transition) or a specific 
heat anomaly (super-critical behavior) . As with the free 
energy and entropy it is impossible to come to a definite 
conclusion regarding the presence and nature of a phase 
transition from this plot. Of particular importance is the 
specific heat for the van der Waals fluid which is seen to 
be nearly constant and equal to 3/2 as it should |jO"|| . 
There is a small slope for the specific heat over the tem- 
perature range in question: for T = T c /2 Cy = 1.507 and 
for T = 2T C C v = 1-502. This illustrates the technique 
employed here, following equations (|o|) and ((2^), yields 
results that are no better than 0.25% for the quantities 
determined in this paper. 

Finally, Figure 13 shows the constant volume caloric 
curves of T as a function of E for each system at the same 
isochores discussed above. Also shown are the values of 
the temperature and energy of the system, E = F + TS, 
along the boundary of the coexistence region. Just as 
with the behavior of the entropy isochores, the behavior 
of the caloric curves for each isochore shows behavior that 
is consistent with either a first order or continuous phase 
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FIG. 12. The specific heat at constant volume versus tem- 
perature for (a) VI, (b) V2, (c) V3 and (d) the van der Waals 
fluid. Full circles show the boundary of the coexistence re- 
gion. Dotted curves are sub-critical. Solid curves are critical. 
Dashed curves are super-critical. 

transition in a small system. Each isochore shows similar 
behavior, a steep rise followed by a region of shallower in- 
cline followed by a portion which approaches E/A = |T. 
The lack of a flat region, or back bend, in the caloric curve 
is not due to the small size of the system but rather due 
to the system being held at a constant volume. Only for 
isobars will flat regions or back bends be observed in the 
canonical nuclear model. See following section. Again 
the van der Waals system shows very different behavior, 
a steady rise in the temperature as a function of energy. 
And again a three dimensional plot is needed to clearly 
understand the nature of each isochore. 

C. Isobars 

The pressure of the system is held constant and the be- 
havior of various quantities with respect to the system's 
temperature is explored. Here some care should be taken 
with the interpretation of the results. The analysis out- 
lined above is still followed. However, since the pressure 
is an extensive quantity plots such as free energy versus 
temperature are now plots of (F) versus T. Where (F) 
is the mid point in the AF range over which a difference 
such as eq. (^3|) is taken. In the limit of vanishing in- 
terval size, this approximation is accurate. Also because 
the pressure is an extensive variable, it was necessary to 




E = F+TS E = F+TS 

FIG. 13. Isochoric caloric curves for (a) VI, (b) V2, (c) V3 
and (d) the van der Waals fluid. Full circles show the bound- 
ary of the coexistence region. Dotted curves are sub-critical. 
Solid curves are critical. Dashed curves are super-critical. 

allow some small variation in P in order to make a plot 
such as (F) against T. The variation in P was usually 
less than one percent; i.e. AP/P < 0.01. Small changes 
in the amount of variation of P had little effect on the 
analysis presented here. Large changes in the variation 
in P wash out the behavior observed below. 

The isobaric free energy as a function of temperature 
is shown for all systems in Figure 14 as well as the values 
of the isobaric free energy along the boundaries of the co- 
existence region. In all systems there is a back bend in 
the free energy curve for sub-critical isobars. The sub- 
critical isobar also traverses the coexistence region. It 
would be possible to perform the Maxwell construction 
procedure and deduce the critical point from these plots. 
The van der Waals fluid system shows that the critical 
point determined in the construction of the P-V coex- 
istence curves agrees with considerations of isobaric F . 
See Figure 14d. The critical isobar shows a vertical slope 
tangent to the coexistence curve and no back bend. The 
super-critical isobar does not traverse the coexistence re- 
gion and shows no back bend. 

Figure 15 shows the isobaric temperature as a function 
of reduced density. Temperature is plotted as a function 
of reduced density in the spirit of the Guggenheim plot 
filf which shows the universal behavior of several fluids 
near their critical point. On the scales shown here no 
universal behavior is observed. It may be that very near 
the critical point, the coexistence curves for each system 
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FIG. 14. Free energy isobars as a function of the temper- 
ature for (a) VI, (b) V2, (c) V3 and (d) the van der Waals 
fluid. Full circles show the boundary of the coexistence re- 
gion. Dotted curves are sub-critical. Solid curves are critical. 
Dashed curves are super-critical. 




0.6 0.7 0.8 0.9 1 



P/P< 



FIG. 15. Isobars of temperature versus reduced density for 
(a) VI, (b) V2, (c) V3 and (d) the van der Waals fluid. Open 
squares show the spinodal. Full circles show the boundary of 
the coexistence region. Dotted curves are sub-critical. Solid 
curves are critical. Dashed curves are super-critical. 



are identical. The T-p isobars for all systems show the 
expected behaviors: a sub-critical back bending curve 
that traverses the coexistence region giving way to the 
critical isobar, a critical curve with a flat section which 
intercepts the coexistence region at the critical point and 
finally a super-critical isobar which avoids the coexistence 
region altogether. 

The behavior of the isobaric entropy as a function of 
temperature is also just as expected. See Figure 16. 

Figure 17 shows the isobaric caloric curves for each sys- 
tem. As the energy is constructed from the free energy 
and entropy of the system the back bending observed 
in the sub-critical isobars is expected. Also shown is 
the value of temperature and energy along the boundary 
of the coexistence region. Note the difference between 
the isochoric caloric curves shown in Figure 13 and the 
isobaric caloric curves shown in Figure 17. When the 
pressure is held constant all of the canonical systems dis- 
cussed here, including the van der Waals fluid, show sub- 
critical caloric curves with a back bend. No back bending 
is present when the volume is held constant. 

Finally Figure 18 shows the constant pressure specific 
heat, Cp, as calculated from: 

taking as input the isobars shown in Figure 17. The re- 



sults for Cp of these canonical calculations are similar 
to the behavior reported in micro canonical models for 
various systems g], p8|] , [ |32]] , J33|] . The sub-critical iso- 
bars show the remnants of poles with Cp < values in 
between. The critical isobar shows the remnants of a di- 
vergence, and the super-critical isobar shows some peak- 
ing behavior. The lack of poles and divergences is not 
due to the finite size of the systems VI, V2 and V3, but 
rather due to the computational nature of these calcula- 
tions. The calculations for the van der Waals fluid are, 
in effect, for a truly thermodynamic system and the van 
der Waals critical isobar of Cp still shows no true diver- 
gence. The critical isobar of Cp for V2 shows a negative 
value which is due to the computational nature of the 
calculation and the manner in which Cp was calculated 
from E and T values. 

D. Iso-nothing: Variable (P, V, T) 

The nuclear model presented here assumes a system 
enclosed in some volume. An actual excited nucleus is 
not enclosed in a volume. This has the effect of forcing 
a path through the thermodynamic phase space which 
is considerably different from any of the paths investi- 
gated thus far. To bridge the gap between reality and 
tractable calculations, an energy dependent free volume 
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FIG. 16. Isobars of the entropy versus temperature for (a) 
VI, (b) V2, (c) V3 and (d) the van der Waals fluid. Full circles 
show the boundary of the coexistence region. Dotted curves 
are sub-critical. Solid curves are critical. Dashed curves are 
super-critical. 



is assumed in models such as the ones presented in this 
work |]. At low energies the free volume of the system 
is assumed to be nearly constant and vanishingly small. 
At a given energy the system is allowed to expand and 
the free volume increases from near zero. This has the ef- 
fect of tracing a path through the thermodynamic phase 
space of the system off any of the trivial paths along one 
of the axes investigated above. It is possible to examine 
the effects of such a parameterization of the free volume 
as a function of energy with the calculations made here. 
Calculations in this section were performed only for VI. 

Following the same ideas presented in the discussion of 
isobars, the values of F, S and T are determined along 
a path in Vf as a function of E. As Vf changes values 
of F, S and E are picked from the appropriate isochore. 
For example, instead of traveling along a path parallel to 
one of the axes in Figure 2a, e.g. an isochore, isotherm 
or isobar, an energy dependent free volume was chosen so 
that the system evolved through thermodynamic phase 
space on a non-trivial trajectory. See Figure 19. In Fig- 
ure 19a small points show a small sample of the set of 
calculations for F(T, V). Larger points show the val- 
ues of F(V,T) selected for a Vf(E) trajectory described 
above. Figures 19b and c show the same for S(T, V) and 
E(T,V). 

For the purposes of the present analysis, values of 
S(T, V) were used from different isochores: as Vf changed 
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FIG. 17. Isobaric caloric curves for (a) VI, (b) V2, (c) V3 
and (d) the van der Waals fluid. Full circles show the bound- 
ary of the coexistence region. Dotted curves are sub-critical. 
Solid curves are critical. Dashed curves are super-critical. 



values of S(T, V) were selected from the appropriate iso- 
chore. If the change from isochore to isochore is small, 
AVf ~ 0, then the procedure is a good approximation. 
The intervals in Vf for the calculation of F(T, V) for 
this analysis were on the order of two percent of the 
free volume of the system; e.g. when Vp — 500fm 3 , 
AV f ~ 10fm 3 and when Vf = 5fm 3 , AVf ~ O.lfm 3 . 

Figure 20a shows three different paths through ther- 
modynamic phase space: one which travels very near to 
the critical point and into the coexistence region (solid 
curve), and two others which avoid the coexistence region 
all together (dotted and dashed curves). As mentioned 
previously for the full version of this model jjj the free 
volume is nearly zero for the lower end of the energy 
range. At some energy, 2 MeV/A for the solid curve in 
Figure 20, the system is allowed to expand and the free 
volume increases as a function of energy. The functional 
form of Vf(E) is not identical to other models W but 
close enough to show the same behavior observed in the 
full version of the model [[!). Also shown in Figure 20a 
are the values of the free volume and energy along the 
boundary of the coexistence region. Note that for the 
solid curve the Vf(E) trajectory enters the coexistence 
region near the critical point and leaves the coexistence 
region at a higher energy and free volume. The other two 
trajectories shown in Figure 20 have the same general 
behavior: increasing free volume with increasing energy, 
but the precise paths differ. 
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FIG. 18. Specific heat at a constant pressure versus en- 
ergy for (a) VI, (b) V2, (c) V3 and (d) the van der Waals 
fluid. Dotted curves are sub-critical. Solid curves are criti- 
cal. Dashed curves are super-critical. 

In a plot of free volume against temperature back bend- 
ing is observed for two of the trajectories presented here. 
See Figure 20b. The solid curve trajectory of Vf(T) be- 
gins with a small free volume that is constant until a 
temperature of just over 8 MeV, then the system expands 
and cools. This is shown by the back bend. The solid 
curve Vf(T) trajectory then enters the coexistence region 
near the critical point. At a temperature between 7 and 
7.5 MeV the slope of the free volume nearly diverges and 
then changes in sign. After this point, further expansion 
in the free volume is accompanied by an increase in the 
temperature of the system. The Vf (T) curve then leaves 
the coexistence region at T < T c and V > V c . The other 
trajectories show similar, but less extreme behavior. 

Other projections of the phase diagram for this model 
with the solid curve Vf(E) trajectory show back bends 
as well. See Figures 20c and d. The Vf-P projections 
shows the system's pressure increases nearly an order of 
magnitude over the constant free volume section. When 
the system is allowed to expand, the pressure drops. The 
trajectory passes near the critical point as it enters the 
coexistence region and in the course of back bending exits 
the coexistence region. 

The solid curve P-T trajectory is equally interesting. 
As the system increases in temperature with a fixed free 
volume, the pressure increases. When the system reaches 
a temperature between 8 and 8.5 MeV the expansion sets 
in and the pressure and temperature both drop so that 




FIG. 19. Solid circles show the path of an expanding nu- 
clear system, VI, through thermodynamic phase space. Nei- 
ther volume nor pressure is held fixed. Points show the cal- 
culations of the free energy surface as shown in Figure 3. 

the trajectory moves towards the critical point. At a tem- 
perature between 7 and 7.5 MeV, the system reverses the 
trend and both pressure and temperature increase. See 
Figure 20d. Again the other two trajectories show similar 
behavior, but to a lesser extent as the Vf(E) trajectory 
becomes smoother. 

As back bending has already been observed for other 
Vf(E) trajectories, it is no surprise that the caloric curve 
for this changing volume system also shows a back bend. 
See Figure 20e. The solid caloric curve shown here is rem- 
iniscent of other back bendin g ca loric curves already pub- 
lished in reference JjJ , [^0| - |22|] where variable free vol- 
ume constrained canonical and constrained grand canon- 
ical calculations are made, but not those published in 
reference Q , [ p8[ where a constant volume micro canoni- 
cal calculation is made. Also shown in Figure 20e are the 
values of the temperature and energy along the coexis- 
tence curve. Again the trajectory of the solid curve vari- 
able free volume enters the coexistence region near the 
critical point and exits the coexistence region at T < Tc 
and E > E c . The dotted caloric curve also shows back 
bending, albeit to a more limited extent and the dashed 
caloric curves shows no major back bend. 

Next the specific heat of the system was calculated via 
eq. (|26|). See Figure 20f. The application of eq. ( pr^ ) 
to this trajectory through thermodynamic phase space is 
problematic. From Figure 20d it is clear that pressure 
is not a constant and thus eq. (|2^) should not be used. 
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FIG. 20. The path of an expanding nuclear system, VI, 
through thermodynamic phase space: (a) the free volume as 
a function of energy, (b) free volume as function of tempera- 
ture, (c) free volume as a function of pressure, (d) pressure as 
function of temperature, (e) the resulting caloric curve, (f) a 
specific heat based on the derivative of the caloric curve and 
(g) a specific heat along the phase space trajectory. Full cir- 
cles show the boundary for the coexistence region. Different 
line types show different (P, V, T) trajectories. See text for 
details. 



However, it has become commonplace to follow this sort 
of procedure [^0| , [p2l even though it is in contradiction 
with definition of Cp or Cy ■ When eq. ( |26| ) is applied to 
the solid and dotted curves in Figure 20e, the resulting 
specific heat shows negative values and the remnants of 
a divergence. 

Finally the specific heat of the system was determined 
in the same manner that the entropy of the system was 
determined. The specific heat at a constant volume was 
calculated along an isochore, then the values of Cy were 
selected from the paths through thermodynamic phase 
space of the V/(E) trajectory. See Figure 20g. For the 
solid like path, the value of Cy shows a steady rise until 
an energy of 2 MeV/ A and then a sharp rise as the system 
expands. As the system continues to expand and the 
energy increases the value of Cy reaches a maximum and 
then shows a gradual decline. No Cy < is observed in 
this plot. The other two paths show smoother behavior. 

The question that now arises is what, if any, insight 
into the nature of the phase transition can be obtained 
from curves such as those in Figure 20e and f. Were 
there no other knowledge of the system, the back bends 



observed in the solid and dotted caloric curves would sug- 
gest that the system had gone through a first order phase 
transition. Negative values and a peak in the specific 
heat would seem to confirm this. While for the dashed 
curve, the lack of back bending in the caloric curve and 
the lack of a negative specific hear would argue either 
for a continuous phase transition of no phase transition. 
However, from the analysis of the previous sections, the 
location of the critical point and the shape and location 
of the coexistence curve are known. The addition of this 
knowledge makes it clear that the naive analysis of Fig- 
ures 20e and f can provide misleading results. While the 
solid caloric curve shows a back bend, the trajectory of 
the systems goes through the critical point, into the co- 
existence region and exits along a sub-critical path. Can 
one conclude that the system has undergone a contin- 
uous phase transition, or has it undergone a first order 
phase transition because the trajectory traverses the co- 
existence region? The answer to the first question is yes, 
since the system does reach the critical temperature and 
density simultaneously The answer to the second ques- 
tion also appears to be "yes" since the specific heat found 
from eq. ( p6| ) is less than zero (for the solid curve) . Fur- 
thermore, the solid curve intersects the coexistence curve 
at a T less than T c . Note, however, that without knowl- 
edge of the location of the coexistence curve, the above 
questions cannot be unambiguously answered. A naive 
inspection of the dotted caloric curve may lead one to 
conclude that the back bending is indicative of a first 
order phase transition. Such is not the case, Figure 20e 
shows that the dashed line never traverses the coexis- 
tence region. It is clear that drawing conclusions based 
on caloric curves is difficult unless one has knowledge of 
the complete thermodynamics of the system. 



E. Critical exponents from thermodynamic 
quantities 

Any model that attempts to describe a system capa- 
ble of undergoing a continuous phase transition should 
exhibit quantities with singular behavior that, near the 
critical point, are described by power laws with a consis- 
tent set of critical exponents. These critical exponents 
should obey well known scaling laws and may, or may 
not, fall into one of the established universality classes. 
To that end, four critical exponent values arc determined 
and three scaling laws are checked for these models using 
the critical point, (T c , P c , p c ), determined previously and 
other thermodynamic quantities. Note that in the deter- 
mination of critical exponents presented here, thermody- 
namic variables are used explicitly, e.g. in the extraction 
of 7 it is the isothermal compressibility that is used, and 
not moments of the fragment distribution. 
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1. Power law results 

The exponent a is determined by the behavior of the 
specific heat along the critical isochore, see Figure 12. 
The CV(T) curve was fit with the functional form: 



Cv(T) = H± 



T-T c 



1 



V 



a± 



G± (27) 



The fit 



on both sides of the critical point T^T C [| 
parameters H± , a± and G± were allowed to vary to min- 
imize the x 2 °f t ne fit- Figure 21 shows the results for 
each system and Table II lists the extracted exponents. 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 



U3.5 
3 
2.5 
2 
1.5 
I 

0.5 







L (a) I 




f / 

r / 

r / 









O 7r 




0.2 0.4 



0.2 0.4 
(T-TJ/T c 





V 


~~~~ 





0.2 0.4 
(T-T c )/T c 



-1 ' 1 1 ' 1 ' 1 ' 1 
: (d) 




j 


L 



-0.5 -0.25 0.25 0.5 

(t-t yr 



FIG. 21. Determination of the critical exponent a from the 
specific heat at a constant density, p = p c for (a) VI, (b) V2, 
(c) V3 and (d) the van der Waals fluid. Solid lines show a 
sample fit. 

The functional form in eq. ( p7|) did not fit the curves 
shown in Figure 21a, b and c over the entire range of 
(T — T c )/T c . To some degree this is to be expected. 
Near the critical point the finite size effects, which mani- 
fest themselves first in the smoothing of the kink of sub- 
critical free energy isochores, diminish a diverging specific 
heat into a peaking specific heat. Far from the critical 
point, the analytic terms in the expression for the spe- 
cific heat become dominant and the power law behavior 
is overwhelmed. In some mid-range region, neither toO 
far from nor too near to the critical point behavior con- 
sistent with eq. (|27j) was observed. Various fits were 
tried on both sides of the critical point, but only those 



which gave a matching value for a were considered. Fig- 
ures 21a, b and c show the results of one such fit. Table 
II lists the average results for many such fits. 

The van der Waals fluid shows much different behavior 
that do the nuclear models. The constant value of Cy = 
3/2 in the van der Waals fluid leads to the result of a = 
as expected for a mean field model. Based on this result 
it would seem that the nuclear models are not mean field 
models. They show a peaking in the specific heat that is 
inconsistent with the behavior of a van der Waals fluid 
type of mean field model or the behavior of the Landau 
model which shows a discontinuity in the specific heat. 

The exponent (3 is determined using the (P, VA points 
along the coexistence curve, shown in Figure 15, which 
should be described by 

P 

(28) 



Pi ~ Pg 



T r -T 



Tr 



Fitting pi — p g versus (T — T c ) /T c to a simple power law 
for the positive slope portion of Figure 22a, b and c gives 
the exponent /3. See Table II for results. The van der 
Waals fluid recovers the mean field value of /3 — 1/2. The 
model VI gives the least impressive fit results. 




10 
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FIG. 22. Determination of the critical exponent /3 from the 
liquid-gas density difference along the coexistence curve for 
(a) VI, (b) V2, (c) V3 and (d) the van der Waals fluid. Solid 
lines show a sample fit. 

Near the critical point the isothermal compressibility, 
kt, is given by 

rp _ rp 7± 

^T = r± (29) 



1G 



TABLE II. Critical exponents for the models 



System 


a± 


13 


7 


5+ 


5+ 


T 


VI 


0.6±0.2 


0.5±0.2 


1.0±0.2 


1.2±0.3 


1.4±0.3 


2.49 ± 0.03 


V2 


0.9±0.2 


0.6±0.1 


1.0±0.4 


2.8±0.1 


2.7±0.2 


2.10 ±0.02 


V3 


1.0±0.2 


0.8±0.1 


1.2±0.2 


2.5±0.2 


3.2±0.2 


2.22 ±0.02 


vdW 


0.0 


0.5±0.1 


0.98±0.02 


2.7±0.2 


3.1±0.3 


2.33 ±0.02 



For T < T c fitting k t versus \(T-T C )/T C \ along the 
coexistence curve gives 7_ while 7+ is determined by 
fitting k t versus |(T - T c )/T c \ for T > T c at p = p c . Due 
to the imprecise nature of the data along the coexistence 
curve for T < T c fits were made only for T > T c . Fits 
were made over the entire region of |(T — T c )/T c \ of T > 
T c . The results for the extraction of the exponent 7 are 
shown in Figure 23. 






KT-rj/Tj 

FIG. 23. Determination of the critical exponent 7 from the 
isothermal compressibility, kt, for (a) VI, (b) V2, (c) V3 and 
(d) the van der Waals fluid. Solid circles show the behavior 
of kt for T > T c and at p = p c . Solid lines show a sample fit. 
Open squares show the behavior of kt along one edge of the 
coexistence boundary. Dashed lines show the fit for T > T c 
multiplied by a constant. 

For VI two different power law regions appear to be 
present, one close to the critical point and one further 
from the critical point. See Figure 23a. The error bars 
on the 7-value in Table II account for this behavior. The 
fit to the entire |(T — T c )/T c \ region is used because the 
resulting power law shows some agreement with the be- 
havior of kt for T < T c when the coefficient of the power 



law is increased by some factor. See dashed line and open 
squares in Figure 23a. Similar arguments apply to the 
results for V2 and V3. See Figures 21b and c. 

The van der Waals fluid shows the expected behavior 
and recovers the value of 7 = 1 to within error bars. See 
Figure 23d and Table II for results. The T < T c behavior 
of kt also shows the expected power law behavior with 
the appropriate exponent value. 

Examining the shape of the critical isotherm leads to 
an estimation of the exponent 8 from: 



\P - P 



<5± 



(30) 



The critical isotherm was examined independently for 
p < pc, which gives S-, and p > p c , which gives <5+. 
As with the exponents a and 7, a system with a continu- 
ous phase transition the values of S± should be the same 
on both sides of the critical density. This fact is again 
used as guide in searching for fitting regions to extract 
the exponent 5. See Figure 24 and Table II. 

For VI only the regions closest to the critical point 
gave matching 5-values. The error bars on the 5-values 
in Table II reflect the changes in 5± when different fit 
regions are examined. In V2 there are regions on both 
sides of the critical point which yield a matching set of S± 
values. No such region could be found for V3, even very 
close to the critical point. The van der Waals fluid shows 
some regions on both sides of the critical point where S± 
match, to within error bars, and agree with the expected 
value of 6 = 3. 

Finally, the topological exponent, r, from Fisher's 
droplet model J35| can be recovered based on considera- 
tions of the compressibility factor, C/ via the relationship 



C 



f 



C(r) 



C(r-l) 



(31) 



The Riemann £ functions of eq. ( [31] ) were summed from 
1 to 1000000000. When a value of r = 7/3 was input 
for the van der Waals fluid, eq. (31) yielded a value 
of 0.393 indicating that terminating the summation at 
1000000000 yields a value of C/ that is approximately 
5% too high; for the van der Waals fluid C f = 3/8. This 
supposition is supported by decreasing the upper sum- 
mation limit and observing and increase in the value of 



17 




ip-pj/p c ip-p t yp< 



FIG. 24. Determination of the critical exponent 5 from the 
critical isotherm for (a) VI, (b) V2, (c) V3 and (d) the van 
der Waals fluid. Open circles show the critical isotherm for 
p > p c and at p — p c . Solid lines show a sample fit. Open 
squares show the critical isotherm for p < p c . Dashed lines 
show a sample fit. 

Cf. This error was accounted for in the estimation of 
the value of r. See Table II for results. 

2. Scaling laws 

With four critical exponents determined it is possible 
to perform a consistency check using the well known scal- 
ing relations. For example, the Rushbrooke inequality 
shows that: 

a + 2/3 + 7 = 2, (32) 

here shown as an equality in keeping with the scaling 
hypothesis and renormalization |f37f . And the Griffiths 
equality: 

a + 13(1 + 6) = 2 (33) 
and the Widom equality: 

/3(«S - 1) - 7 = 0. (34) 
And finally from Fisher's droplet model: 

^-P^ = 0. (35) 
7 3 — r 



Using the average values determined for a, f3, 6, 7 and 
r the results for these scaling laws are compiled in Ta- 
ble III. Only the van der Waals fluid results consistently 
satisfy the above scaling laws to within error bars. The 
nuclear models generally fail to satisfy three of the four 
scaling laws. This failure is inconsistent with the behav- 
ior of the phase diagram, shown for example in Figure 5, 
which appears show a critical point, thus indicating the 
presence of a continuous phase transition. While mod- 
erately good fits are observed for the specific heat, the 
liquid- gas density difference, the isothermal compressibil- 
ity and the critical isotherm for each of the versions of 
the nuclear model, the meaning of these power laws and 
critical exponents remains an open question in light of 
the failure to adhere to well known scaling laws. 

V. SUMMARY 

It has been shown that the type of nuclear model dis- 
cussed here exhibits many features commonly associated 
with a system in which critical phenomena are present, 
e.g. a coexistence curve, power laws, critical expo- 
nents. By removing both the Coulomb and temperature- 
dependent free energy terms, it was found that the ap- 
pearance of a critical point in these models is due to the 
interplay between the surface, volume, and translational 
free energy terms. However, these types of models are 
not without inconsistencies. One striking inconsistency is 
the fact that the temperature dependent surface free en- 
ergy gives rise to a infinite negative specific heat at the 
critical temperature used by the model. Furthermore, 
no version of the model showed a critical temperature 
that agreed with the one explicitly input into the sur- 
face free energy term. Additionally, when the surface 
term was rendered temperature independent the critical 
point remained. Thus suggesting that the appearance of 
a critical point in these models is not dependent on the 
temperature dependence of the surface term but rather is 
a result of the interplay between the surface and volume 
free energy terms. 

The critical temperature and density have been deter- 
mined by examining isotherms in the P — p plane. In 
the neighborhood of this critical point, singular behav- 
ior characterized by power laws was observed. However, 
these critical exponents do not obey well known scaling 
relations. This is a particularly troublesome occurrence 
as any model with true critical behavior, even the sim- 
ple van der Waals fluid, does have exponents which obey 
these scaling relations. It is possible that an examina- 
tion of this model for larger systems, with smaller steps 
in temperature and volume in the calculation of the free 
energy, will yield a consistent set of critical exponents. 

It is important to note that the critical densities found 
here are much higher than could be realized with a clos- 
est packing of normal density nuclei. Additionally, these 
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TABLE III. Scaling law results 



System Rushbrooke Griffiths Widom Fisher 

VI 2.6±0.4 1.8±0.5 0.9±0.2 -0.5 ± 0.2 

V2 3.2±0.5 3.2±0.4 -0.1±0.2 0.5 ± 0.3 

V3 3.7±0.3 3.9±0.6 -0.2±0.3 0.4 ±0.1 

vdW 2.0±0.1 2.0±0.4 -0.0±0.2 0.0 ± 0.1 



critical densities are significantly higher than those typi- 
cally used to compare model predictions 0, to data. 

A major conclusion of this work is that the particular 
phenomenological description of the free energy of a hot 
nucleus leads to several inconsistencies regarding both 
temperature and density. It was pointed out that the 
parameterization of the surface free energy leads to a 
negative and divergent contribution to the specific heat 
at T approaches the value of the parameter T c in eq. 
(|J). Furthermore, all values of the critical temperature 
found from examination of isotherms in the P-p plane 
are much below this parameter value. Thus, while use of 
such a model may well lead to an excellent description of 
multifragmentation data, the lack of internal consistency 
noted here makes the interpretation of data in terms of 
the model problematic. Such agreement may rest more 
on the phase space sampling and variable free volume 
inherent in the model than on the finer details examined 
here. 

Finally, it has been shown that the variable volume 
version of this phenomenological model of multifragmen- 
tation exhibits caloric curves which can be misinterpreted 
in the absence of detailed knowledge of the complete ther- 
modynamic phase diagram. 
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